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ABSTRACT 

Seven cosmological hydrodynamic simulations of disk galaxy formation are analyzed to determine 
the alignment of the disk within the dark matter halo and the internal structure of the halo. We find 
that the orientation of the outer halo, beyond ~ 0.1 r v i r , is unaffected by the presence of the disk. 
In contrast, the inner halo is aligned such that the halo minor axis aligns with the disk axis. The 
relative orientation of these two regions of the halo are uncorrelated. The alignment of the disk and 
inner halo appears to take place simultaneously through their joint evolution. The disconnect between 
these two regions of the halo should be taken into account when modelling tidal streams in the halos 
of disk galaxies and when calculating intrinsic alignments of disk galaxies based on the properties of 
dark matter halos. 

Subject headings: galaxies: formation — galaxies: halos — dark matter — galaxies: spiral — galaxies: 
evolution 



1. INTRODUCTION 

The alignment of galactic disks within their triaxial 
dark matter halos is both observationally and theoreti- 
cally difficult to determine. Observationally, dark matter 
can only be traced by its gravitational effect on lumi- 
nous matter. Efforts to determine the shapes of indi- 
vidual dark matter halos around disk galaxies usually 
rely o n the motions of tidal s tream s llJohnston et al.l 



12003 iThata et al 



Martmcz-Delgado ct al. 



rings (Sackctt ct al 
( Oili ng fc Merrifield 



1994 
200' 



2001, 2004: Hcl mil 12004a 
2001 lLaw et alJ l200l . polar 



or the flaring of the gas layer 
Most of these studies have 



found that the halo is flattened along the pole of the 
galactic disk; however, the degree of flattening for even 
the best-studied case of the Milky Way is controversial, 
with values ranging from 0.8 to 1.7 (i.e. elongated along 
the Galactic pole). 

Theoretical determinations of galactic disk alignment 
within halos have been restricted by available computing 
power. A-body simulations have been used to deter- 
mine the flattening and internal alignment of pure dark 
matter halos l | B arnes fc Efstathiou||1987|: |Warr en et al.l 
Il992t Uing fc Sutoll2002t iBailin fc Steinmetzll2005i here- 
after IBSOM) . These studies have found that the axes of 
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halos, especially the minor axes, are very well aligned 
internally; the halo is well approximated by a set of con- 
centric ellipsoids. Under the assumption that the an- 
gular momentum of the baryons which form the disk is 
aligned with the angular momentum of the dark mat- 
ter, the disk axis is e xpecte d to lie typically 25° from 
the halo minor axis (BS05). However, hydrodynami- 
cal simulations suggest that the angular momentum of 
the baryons and dark mat t er are not perfectly aligne d 
(van den Bosch et"^l200l ISh arma fc Stemmed EPOS). 
Moreover, the orientation of t he halo may change in the 
pres ence of a misaligned disk l)Binnev et al.l lY998. here- 
after HjdH). 

This problem can now be addressed with recent 
high resolution simulations of disk galaxy formation. 
These simulations have traditionally formed disks that 
are much more compact th an observed galactic disks 
({Navarro fc Steinmetzl I1997T) . However, modern treat- 
ments of feedback and more consistent analyses of the 
simulations and observations have led to much better 
agreement: the physical size of the disks is realistically 
reproduced, and the lack of angular momentum mani- 
fests itself primarily in a ve ry pronounced slo wly rotat- 
ing bulge-halo component ijAbadi et al.l 12003^ . Recent 
works employing advanced multi-phase descriptions of 
the ISM reduce the dominance of the bulge, in particu- 
lar if feedback from AGNs is included ({Robertson et alJ 
120041 lOkamoto et al.l200HD . Since the physical size of the 
disks is realistically reproduced and since the location of 
the surrounding matter which provides the tidal torque 
is unaffected, the orientation of the angular momentum, 
a nd therefore th e disk, is expected to be robust. 

INavarro et al.l l)2004f) found that in four such simula- 
tions, the disk axis tends to lie within 30° of the in- 
termediate axis of the large scale mass distribution at 
turnaround, indicating that primordial tidal torques play 
an important role in determining the final disk orienta- 
tion. It is known that the radial distribution of dark 
matter i n the halo changes in the presence of cooled 
baryons (jGnedin et al. 1 12004]) . and therefore it is likely 
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that other p roperties of the matter d istribution a re af- 
fected too. iKazantzidis et al.l (|2004Tl (hereafter IK04I) 
studied the change in the axes of the halo of a high res- 
olution disk galaxy simulation compared to the iV-body 
case. They found that the axis ratios of the halo were 
dramatically higher (i.e. more spherical) when the sim- 
ulation included baryon cooling than when the baryons 
were not present or not allowed to cool to a disk. They 
also found that the minor axis of the halo was aligned 
with the disk axis within r < 0.2 fiso, but that the halo 
became elongated along the disk axis beyond that radius. 

The alignment of disks with their dark matter halos 
is important for understanding a number of unresolved 
questions regarding disk galaxy formation, evolution, 
and dynamics. Warps i n disk galaxies are ubiquitous 
JRe^hetnikpv^^COTnbes! 11993 ISchwarzkopf fc Dettmarl 
l200lITGarc7a^RuTr^r^Lf|2002|) and have been proposed 
to be caused by misalign ment between the disk and 
the dark matter h alo (Dckcl & Shlosman 198 31 iToomrd 
1983HBailirJl2004t see however INelson fc Tremainelll995t 
BJD98). This possibility can be ruled out if such 
misalignments never occur in fully self-consistent sim- 
ulati ons of disk galaxy formation. The Holmberg ef- 
fect llHo lmberd 119741: iZaritsky e t al.lll997l see however 
lBrainerdl2004tlWilhnan et al-lboolT mav be explained if 
galactic disks preferentially lie perpe ndicular to the ma- 
jor axis of their da rk matter halo l)Knebe et alJ 120041 
iZentner et al.l 12005^ . Models of the Sagittarius tidal 
stream cannot simultaneously fit the leading and trail- 
ing arm of the stream with t he same halo flattening 
l|Helmill2004bl lLaw et al.ll2005|) : a more detailed under- 
standing of the relationship between the location of the 
disk and the shape of the dark matter halo may be neces- 
sary to resolve this issue. Finally, the intrinsic alignment 
of observed galaxies acts as a background contaminant 
in weak lensing studies. Predictions of the magnitude 
of this effect based on linear theory or TV-body simula- 
tions must assume a disk orientation for each halo, which 
is as yet untested llHeavens et al.ll2000t iPen et aLll2000t 
ICrittenden et a,lJl2001ft. 

In this Letter, we analyze the shape and internal align- 
ment of the dark matter halos of seven high resolution 
cosmological disk galaxy formation simulations, and dis- 
cuss the possible observational effects of such internal 
structure. 

2. METHODOLOGY 

We analyze seven cosmological simulations which use 
the multi-mass technique to self-consistently model the 
large scale tidal field while simulating the galactic disk 
at high resolution. These simulations include self- 
consistently almost all the important physical processes 
in galaxy formation, such as self-gravity, hydrodynamics, 
radiative cooling, star formation, supcrnovae feedback 
and metal enrichment. 

The properties of the simulations are listed in Table 
Column 1 contains the name by which we refer to each 
simulation, Column 2 is the code used to run the sim- 
ulation, Column 3 is the virial mass of the galaxy, Col- 
umn 4 is the virial radius of the galaxy, and Column 5 
contains the radial extent of the gas disk in each sim- 
ulation, defined as the largest radius at which we find 
gas particles in the disk plane. The primary simula- 
tion we analyze is the "KGCD" simulation, performed 



TABLE 1 
Properties of simulations 



Simulation Name 


Code 


M vir 


^*vir 


»"di s k 






(M ) 


(kpc) 


(kpc) 


KGCD 


GCD+ 


8.8 x 10 11 


240 


10 


AGCD 


GCD+ 


9.3 x 10 11 


270 


21 


AGSPH 


GRAPESPH 


9.4 x 10 11 


270 


23 


KIA5 


GRAPESPH 


1.7 x 10 12 


300 


4 


KIA9 


GRAPESPH 


2.1 x 10 12 


320 


20 


KIB1 


GRAPESPH 


2.9 x 10 12 


390 


16 


KIB2 


GRAPESPH 


8.5 x 10 11 


260 


5 



using the GCD+ code l)Kawata fc Gibson! 12003]) . The 
si mulation is a higher resolution model of galaxy "Dl" 
in lKawata et al.l (|2004). The mass and softening length 
of individual gas (dark matter) particles in the highest- 
resolution region are 9.2 x 10 5 (6.2 x 10 6 ) M and 0.57 
(1.1) kpc, respectively. We analyze the KGCD simula- 
tion both at the final output at z = 0.10 and an output 
at z = 0.37 (approximately 2.2 Gyr earlier) at which the 
disk first appears fully formed in approximately its final 
state to determine whether there is significant evolution 
in the halo structure in the absence of major mergers. 
In order to ensure that our results are not particular to 
this code, the feedback prescription, or the particular 
set of initial conditions, we also an alyze the fina l snap- 
shot of the simulation described in lAbadi et alJ ( 2003) 
("AGSPH"), a simulation with the same initial condi- 
tions as AGSPH but run usin g GCD+ ("AGC D"), and 
four additional GRAPESPH ijsteinmetzl Il996j) simula- 
tions: KIA5, KIA9, KIB1, and KIB2. Dark-matter-only 
versions of each of these simulations were also run in or- 
der to directly compare the alignment of each halo with 
and without a galactic disk. 

The principal axes and axis ratios of t he dar k matter 
halo are calculated using the method of IBS05L Briefly, 
the particles are split into spherical shells bounded by 
radii of r = 0.025,0.05,0.1,0.25,0.5, and 1.0 r vir n . 
The reduced inertia tensor is calculated for the particles 
within each shell and diagonalized to the determine the 
principal axes. Finally, the axis ratios determined from 
the eigenvalues of the reduced inertia tensor are corrected 
for the bias introduced by using spherical shells. In this 
way^ the results at each radius arc entirely independent 

3. RESULTS 

The change in axis ratio for our simulations is shown in 
Figure n An increase in sphericity is clearly evident in 
this figure, with a magnitude that decreases at larger 
radii. The magnitude of this change is smaller than 
that found by K04|, who also found that the presence 
of baryonic cooling dramatically increases the axis ratios 
of dark matter halos, with changes in b/a and c/a that 
range from 0.35 at small radii to zero at the virial ra- 
dius. The change in axis ratios was particularly evident 
in their galactic disk simulation, where b/a rose by almost 

11 We use the fi tting function from Appendix A of 
IKitavam a & Suto ( 1993) to calculate the overdensity that defines 
the virial radius r v ; r as a function of cosmology and z 
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Fig. 1. — Axis ratios b/a and c/a of the dark matter halos 
in simulations including baryonic physics compared to simulations 
containing only dark matter. Blue symbols and lines (dashed liens) 
represent the axes measured at 0.1 r v i r , while red symbols and lines 
(solid lines) represent the axes measured at 0.5 r v ^. Filled symbols 
represent the simulations containing baryonic disks, while open 
symbols represent the ./V-body only simulations. Lines connect 
simulations with the same initial conditions and show the change 
in axis ratios induced by the inclusion of baryonic physics. The 
different symbols refer to the different simulations. 
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Fig. 2. — The alignment between the minor axis of the dark 
matter halo at different radii and the disk axis. The different colors 
and symbols represent different simulations. 



0.6 when a galactic disk was present. This difference in 
magnitude may simply be small number statistics: the 
./V-body version of the galactic disk simulation of K04 
is exceptionally flattened, and therefore a much greater 
increase in axis ratio is possible. 

The direction cosine between the disk axis and the halo 
minor axis at each radius is shown in Figure [21 while Fig- 
ureJSJshows the change in the alignment of the minor axis 
at each radius compared to the dark-matter-only version 
of the same simulation. This is unity if the axes are per- 
fectly aligned and vanishes if they are perpendicular. The 
disk axis is defined to be the minor axis of the gas dis- 
tribution within rdisk! in all cases, this matches the axis 
determined by visual inspection, and that determined in- 
cluding all baryons within rdi s k- Two regimes are clearly 
visible in Figures El and |21 

1. r < 0.1 r V i r : In the inner regions of the halo, the 
minor axis is extremely well aligned with the disk 
axis. 




> 0.1 



In the outer regions of the halo, the 



Fig. 3. — The alignment between the minor axis of the dark 
matter halo of each hydrodynamic simulation and of the matching 
./V-body simulation, as a function of radius within the halo. The 
different colors and symbols represent different simulations. 



orientation of the minor axis is unchanged from the 
iV-body only case. 

The transition between these regimes occurs at a slightly 
larger radius in the AGSPH and KIA5 simulations; as 
these simulations contain respectively the largest and 
smallest disks, the transition radius does not appear to 
depend on the disk radius. 

KIA5 appears as an outlier on all of these plots: it has 
a large change in axis ratio even at large radius, it shows 
very little change in orientation between the inner and 
outer halo, and the minor axis in the hydrodynamic sim- 
ulation is not at all aligned with the iV-body case. This 
latter point may be because the iV-body version of KIA5 
is exceptionally prolate in its outer regions, making a de- 
termination of the minor axis ill-defined. However, we 
note that KIA5 contains the smallest gaseous disk and 
may represent an earlier-type galaxy than the other sim- 
ulations, with a disk that is much less important to its 
dynamics. This interpretation is bolstered by the fact 
that KIB2, which contains the second-smallest disk af- 
ter KIA5, also shows smaller than average misalignment 
between its inner and outer regions. 

The distribution of cosines between randomly-oriented 
unit vectors is uniform with a mean of 0.5. This is pre- 
cisely the distribution of cosines between the inner and 
outer halo axes in Figure|2 Therefore, we conclude that 
the orientation of the inner and outer halo in simula- 
tions containing galactic disks are uncorrelated. In con- 
trast, the minor axes of pure dark matter halos are very 
well aligned throughout their entire extent, with a direc- 
tion c osine of always greater than 0.88 (see Figure 8b of 
BS05J. Therefore, the presence of the disk strongly mod- 
ifies the shape of the inner halo, reorienting it so that the 
halo minor axis aligns with the disk axis. This is in qual- 
itative agreement with IB.ID98L who discussed the effect 
of misalignment on the disk rather than the halo, but 
concluded that the halo reorients itself to line up with 
the disk axis. The results for the KGCD simulation at 
z = 0.10 and z = 0.37 are very similar, indicating that ei- 
ther this reorientation must happen on a short timescale, 
or that the processes that determine the disk and halo 
axes act on both simultaneously during their formation. 

There is some relation between the disk axis and the 
minor axis of the unperturbed halo, as the distribution of 
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innermost points in Figure[!|]has a mean greater than 0.5; 
this is expected due to the tendency of th e halo angular 
momentum to align with its minor axis l) Warren et al.l 
1992; BS05). However, this correlation is not strong and 
the slight misalignment within pure dark matter halos 
ensures that there is no residual correlation between the 
inner and outer regions of halos containing disks. 

4. DISCUSSION 

The final orientation of the galactic disk in a cosmo- 
logical dark matter halo is still unresolved. It is clearly 
related to the angular m omentum of the mater ial in the 
halo, and the results of iNavarro et alJ l)2004fl indicate 
that the primordial tidal torques of the large scale mat- 
ter distribution play an important role, but there is con- 
siderable galaxy-to-galaxy variation in the results of the 
simulations. While the primordial tidal torques generate 
the angular momentum in the material which eventu- 
ally becomes the disk, the accretion of t his material by 
the d isk is a clumpy, stochastic process fVitvitska et, al.l 
2002) , and therefore we are not able to predict the precise 
orientation of the disk. 

However, the results presented in § |3 clearly indicate 
that the halo orientation within 0.1 r v ir matches the disk 
axis. Beyond 0.1 r V j r , the halo orientation is unaffected 
by the presence o f the disk, although the axis ratios are 
somewhat larger l)K04fl . This outer region is itself inter- 
nally well-aligned, resulting in two distinct regions of the 
halo with unique and uncorrelated orientations. These 
results are independent of the time that the halo is stud- 
ied, the code used, and the particular galaxy simulated. 

It is interesting to consider whether the disk drives the 
orientation of the inner halo or vice versa. While the disk 
is less massive than the inner halo (the baryons represent 
15 - 42% of the mass within 0.1 r v ; r ), it is also much 
more flattened than the halo. To evaluate which is most 
important, we have constructed rings of test particles 
tilted by 10° with respect to the disk plane and compared 
the gravitational torque on these test particles due to the 
dark matter with that due to the baryons. For rings of 
radius less than ~ 0.1 r V j r , i.e. in the aligned region, 
the torques from the dark matter and baryons are of 
comparable magnitude, while the dark halo dominates 
the torque at larger radii. Therefore, it appears that the 
alignment is due to the simultaneous evolution of the disk 
and halo rather than the disk directly driving the halo 
orientation or vice versa. 

If the structure of the stellar halo is determined by the 
dark halo, we may expect the properties of the stellar 
halo to change abruptly at 0.1 r v ; r . In fact, the flattening 
of the stellar halo of the Milky Way changes from a flat- 
tened distribution aligned with the disk at R < 15 kpc 
to an essentially sphe rical distribution at R xs 20 kpc 
l)Chiba fc BeerslEOQOD . This spherical distribution may 
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be an intermediate stage between two misaligned flat- 
tened regions. However, the virial radius of the Milky 
Way is likely at least 250 kpc, indicating that this ob- 
served transition occurs at a smaller radius than the tran- 
sition seen in our simulated halos. 

Because the disk is always aligned with the halo in the 
inner regions, general misalignment between the disk and 
halo is not a viable mechanism for generating galactic 
warps. However, if very low surface density disk mate- 
rial existed beyond 0.1 r v ; r , it may settle into the outer 
halo symmetry plane, and theref ore be w arped with re- 
spect to the inner disk. Indeed, IB.ID98I suggested that 
the warp modes within live halos would be essentially 
different from those in static halos not in the configura- 
tion of the disk, but rather in the configuration of the 
halo. However, the outer extent of the baryonic disk 
is determined by the material which has most recently 
cooled; this material is esse ntially unaffected by the an - 
gular momentum problem ijNavarro fc Steinmetd h.997'1. 
Therefore, the maximum extent of simulated disks is ex- 
pected to be accurate, and therefore the disk is unlikely 
to extend into the misaligned outer halo. However, sim- 
ulations of such systems should be carried out to deter- 
mine if such warped configurations are stable, and to 
determine how the surface density of the disk affects the 
tendency of the halo to realign itself. 

The solut i on to the Holmberg effect proposed by 
iKnebe et alJ l)2004|) requires that galactic disks are pref- 
erentially perpendicular to the major axis of the halo. 
However, we find that the orientation of the disk is uncor- 
related with the orientation of the major axis of the outer 
halo, where most satellite galaxies lie, casting doubt on 
this hypothesis. 

Finally, we note that the orbit of the Sgr dSph in the 
Milky Way con tains segments both in side and outside 
of 0.1 r vir (e.g. iHelmi fc White! 12001ft . If the orienta- 
tion of the outer Milky Way halo is significantly different 
from the orientation of the Galactic disk, as suggested by 
our results, then the orbits of particles in the Sagittarius 
stream may differ from those calculated assuming that 
the halo can be approximated by concentric ellipsoids 
(note, however, that the stream depends on the flattening 
and orientation of isopotential contours, which are gen- 
erated by matter at a range of radii and therefore change 
more smoothly than the density). This additional ingre- 
dient may be required to self-consistently model both the 
leading and trailing arms of the stream. 
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